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Abstract 

A general program to show quantum-classical correspondence for bound con- 
servative integrable and chaotic systems is described. The method is applied 
to integrable systems and the nature of the approach to the classical limit, 
the cancellation of essential singularities, is demonstrated. The application to 
chaotic systems requires an understanding of classical Liouville eigenfunctions 
and a Liouville spectral decomposition, developed herein. General approaches 
to the construction of these Liouville eigenfunctions and classical spectral pro- 
jectors in quantum and classical mechanics are discussed and are employed to 
construct Liouville eigenfunctions for classically chaotic systems. Correspon- 
dence for systems whose classical analogs are chaotic is discussed, based on 
this decomposition, in the following paper 0. 

I. INTRODUCTION 

The correspondence principle requires the laws of quantum dynamics to reproduce the 
predictions of classical mechanics when the energy and mass of a dynamical system are very 
large with respect to Planck's constant ||. In other words the correspondence principle 
asserts that the fundamental laws of physics which are valid on the atomic scale are also 
responsible for the observed classical dynamics of macroscopic systems. It is therefore sur- 
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prising, considering the important conceptual role of the correspondence principle that there 
is little direct evidence to support its general validity |3|-§]. In fact it has been controver- 
sially |5]|| argued that the absence of chaos in quantum dynamics plus its presence in the 
dynamics of classical systems prevents any possibility of correspondence || . The main goal 
of this paper and its companion |I| is to provide a consistent theory of correspondence in 
both integrable and chaotic systems. 

To understand correspondence requires that we first adopt a framework in which quan- 
tum and classical dynamics can be realistically compared. The Liouville picture in classical 
mechanics in combination with a phase space representation of quantum dynamics provides 
such a framework. Here both classical and quantum mechanics are represented by objects 
defined on phase space and by operators defined on these objects 0. Our correspondence 
approach is then based on the desire to prove that the quantum dynamics of phase space 
distributions obeys classical mechanics in the classical limit. We do so by focusing on the 
correspondence of what we argue to be the essential elements in dynamics, the eigendistri- 
butions and eigenvalues of the classical and quantum Liouville operators [^,0 (the Poisson 
bracket with the classical Hamiltonian H, and the commutator with the Hamiltonian oper- 
ator H, respectively). 

We outline this correspondence program in Section |T| and describe the progress made 
and developments required in the application to regular and chaotic systems. The approach 
emphasizes the significance of the eigenfunctions of the classical Liouville operator and of 
the classical spectral projectors. These eigenfunctions are obtained in Section |T| for both 
regular and chaotic systems. The correspondence program is then implemented for integrable 
systems in Section [IV| and the nonsingular nature of the approach to the classical limit is 
displayed. This approach to correspondence proves inappropriate for chaotic systems. In 
this case a description of the dynamics in terms of classical spectral projectors is required. 
These quantities are introduced and discussed in Section [V|. A companion paper U shows 
classical-quantum correspondence for chaotic systems. 
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II. LIOUVILLE PICTURE: THE CORRESPONDENCE PROGRAM 



Arguments in favor of the Liouville picture as the conceptual framework most suited for 
the study of correspondence have been presented elsewhere flOlJlll . Here we briefly summa- 
rize. First, the Liouville picture of classical dynamics, which deals with the time evolution 
of phase space densities, is conceptually superior to the trajectory viewpoint because it im- 
plicitly recognizes the imperfect nature of classical preparation and measuring devices. That 
is, one can not initially prepare the exact points in phase space which define a single trajec- 
tory. Similarly, the limitations of quantum mechanical preparation and measurement can be 
incorporated in the von Neumann (quantum Liouville) equation for the density matrix, i.e., 
the von Neumann equation admits solutions which are not states of maximal information. 
Thus, the Liouville equation plays a similar role in quantum and classical dynamics. Second, 
our approach takes advantage of the similarity between the dynamics of classical distribu- 
tions in phase space, as governed by the classical Liouville equation, and the dynamics of 
quantum density matrices, governed by the von Neumann equation in the Wigner-Weyl rep- 
resentation. Adopting this perspective, both classical and quantum mechanics have similar 
structures and rely upon similar quantities. Finally, formal classifications of classical systems 
as integrable, ergodic, mixing, etc. are most properly done in terms of the eigenfunctions 
and eigenvalues of the classical Liouville Operator J[^], i.e., the Poisson bracket with the 
Hamiltonian. Hence this provides the most fundamental of frameworks. 

The time evolution of the quantum density matrix p is given by the solution to the von 
Neumann (or quantum Liouville) equation: 

?l = -i%- 1 [H,p] = -iLp, (1) 

where L = TT l \H , •] is the quantum Liouville operator, and [•, •] is the commutator. Solutions 
to Eq. (0) are of the form: 

P = Px, a e- lXt (2) 
if Px,a satisfies the eigenvalue problem 



Lf>\, a = \p\, a . (3) 

The label a is introduced to accommodate degeneracies associated with the eigenvalue A. In 
particular, a complete characterization of the quantum Liouville eigenvalue problem requires 
the specification of a complete set of commuting superoperators in the same way that solution 
of the Schrodinger eigenvalue problem requires a complete set of commuting observables || . 
This complete set of commuting superoperators can be constructed in the following way. 
Suppose 

H, Ki, . . . , i^r'-l (4) 

for r' < s is a complete set of observables for a system of s degrees of freedom such that 
[H, Ki] = and [Ki,Kj] = for all i,j < r' — 1. It follows that the complete set of 2r' 
commuting quantum superoperators is 

l, n, ^[k lt •], •]+, . . . , \[k r '-i, ■], \[k r ^ •]+ (5) 

where TL = h[H, ■}+ is the energy superoperator, and [•, •]+ denotes the anticommutator, i.e., 
[A, B\+ = Ab + BA. A Liouville eigenstate for the system is completely specified if it is an 
eigenstate of all 2r' superoperators. Thus pa, q in Eq. (|]) is an eigenfunction of the collection 
of superoperators in Eq. (H), with the label a denoting the eigenvalues of all superoperators 
other than L. 

Equations (||) and ([3]) are of the following form in the Wigner-Weyl representation [13 
[x = (p, q) denotes the collection of momenta p and coordinates q] 

dp w (x,t) 
dt 

and 



-zL(xK(x,t) (6) 



L(xK iQ (x) = ApX«(x), (7) 

where 
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2 2 
L(x) = -i#(x) sin(fta/2) = -iHCx) sin 



h d_d_ _ d_d_ 
2 \dq dp dp dq 



(8) 



Here a implies the Poisson bracket, i.e., 



A(x)a£(x) = {A,B} 



&4(x) &B(x) M(x) <9.B(x) 



(9) 



<9q dp dp <9q 

and the arrows over the derivatives indicate that the derivative is taken of the function 
preceding or following the derivative operator. 

Given the solution to Eq. (|7]), the dynamics of any distribution p w (x, t), with initial 
condition p w (x, t = 0) = p(x, 0), can be written as: 



p-(x,t) = £c A , Q p^(x)e 

A, a 



-iXt 



(10) 



where 



dx p(x ,o)p^;(x ). 



(ii) 



That is, 



p w (x,t) = J dx p(x ,0) 



A,« 



-jAt 



(12) 



Consider now the analogous approach in classical mechanics which considers the dynam- 
ics of phase space distributions p(x, t). In particular, p(x, t) satisfies the Liouville equation 

<9p(x, t) 



at 



-iL c p(x, t) 



(13) 



where the classical Liouville operator L c is given by 

L c = z{#(x),-}. 

Here -ff(x) is the Hamiltonian and {•, •} denotes the Poisson bracket [Eq. 
evolution of p(x, t) is given by 



p(x,t) = p A c >a (x)e" 



«A c t 



(14) 

. The time 
(15) 



if Pa c , q ( x ) satisfies the eigenvalue problem 



L c p X c ia (x) = A c pA=, a (x). (16) 

Once again a contains the labeling associated with states which are A c degenerate. In 
particular, a complete characterization of the classical Liouville eigenvalue problem requires 
the specification of a complete set of commuting operators on the classical Hilbert space 
fllljlljl . Suppose the classical analog of our quantum system has r independent constants 
of the motion. Then we can construct r — 1 functions Kj(x) which are constants of the 
motion, i.e., {if (x), Kj(x.)} = 0, and which are in involution so that {^(x), Kj(x)} = for 
all i,j < r — 1. The complete set of commuting operators on the classical Hilbert space is 
then 

L c , if(x), i{#i(x), •}, Ki(x), . . . , i{tf r _i(x), •}, tf r _i(x), (17) 

where L c , i{Ki(x), •},..., i{if r _x(x), •} are first order linear differential operators and 
H(x), Ki(x), . . . , if r _i(x) are multiplicative operators. A Liouville eigenstate for the classi- 
cal problem is usually completely specified if it is an eigenstate of all 2r classical operators 

Given such solutions, the evolution of any classical distribution p(x, t) can be written as 

p(x, t) = J2 c\°,a PA^(x)e- iAC * (18) 

A c ,a 

where 

c\o, a = J dx p(x o ,0)p* A c a (xo), (19) 

or 

(20) 

Both Eqs. ( |1"2"D and ( |20| ) express the time evolution of the phase space density in terms 
of the initial distribution p(x , 0), the eigenvalues A, A c and a term comprised of products of 
Liouville eigenfunctions contained within brackets. This rewriting of the right hand sides of 
the equations, and the term in brackets in particular, will prove to play a central role in the 
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p(x,t) = J dx p(x ,0) 



^p Ac , Q (x)pV iQ (x ) 



X c ,a 



correspondence program we develop in these papers. Indeed the term in brackets, in Eqs. 
(|T2|) and (|20|), is related to the Liouville spectral projection operators, which are heavily 
emphasized in this paper. 

A comparison of Eqs. (|T2| ) and ( p0|) motivates our correspondence approach. Specifically, 



it suggests that one can show quantum-classical correspondence by showing that the h — > 
limit of the quantum Liouville eigenfunctions and its eigenvalues are the classical Liouville 
eigenfunctions and eigenvalues. This approach has proven successful for integrable systems 
10| . However, as shown below and in the following paper Q, it is an incorrect approach 



for chaotic systems. In such cases the quantum Liouville eigenstates have no classical limit. 
Rather, one must deal with quantum and classical spectral projection operators, which have 
singularity-free classical limits. These are related to the objects contained within brackets 
in Eq. (0) and ©. 

The status of this correspondence program is readily summarized for integrable systems 
(r — r' — s). Berry |Tj| demonstrated correspondence for the stationary Liouville eigen- 
functions |n)(n| where |n) are eigenstates of the Hamiltonian operator with corresponding 
energies E n and the s integers n label the quantum states of the Schrodinger picture. Specif- 
ically, he showed that 

plM^h-^J C /ve^(q-v/2|n)(n|q + v/2) ^^ 2 pi n ,o(x) (21) 



in the /i — > limit, where I n = (n+/3)fr, the constants (3 are the Maslov indices [|T^,|T^], p™ n is 



the Wigner transform of the stationary Liouville eigenstate |n)(n|, and the classical Liouville 
eigenfunctions pr,k( x ) axe defined below [Eq. fl30|)1. Correspondence for the nonstationary 
Liouville eigenfunctions |n)(m| and eigenvalues A nim = (E n — E m )/h was considered by Jaffe 



and Brumer |T0[. In particular, they examined the Wigner transform Pn m (x) of |n)(m| and 
showed that for n ^ m 

C.W = h- s/2 J dve^iq- v/2|n)(m|q + v/2) - ^ /2 pi n , m ,n- m (x), (22) 

and 



An.m — > At (23) 

as h — > 0, where they surmised that I njm = (I n + I m )/2. However, as we will later show, this 
In,m is an approximation to the exact result which depends intimately on the eigenvalues of 
the complete set of commuting observables [Eq. (HJ)]. These classical Liouville eigenvalues 
Af, are defined below [Eq. fl3l|)1. 

By contrast, the correspondence program for systems whose classical analogs are chaotic 
has been left largely undeveloped. The nature of the classical Liouville eigenf unctions and 
eigenspectrum is not clearly understood JD], and only the correspondence limits of the 



stationary Liouville eigenfunctions have been explored. Berry |TJ[ and Voros [[HJ were able 



to show that the Wigner function associated with a stationary Liouville eigenstate |n)(n| 
(here r = r' = 1) concentrates in the classical limit to a uniform distribution over the 
classical energy surface of the same energy E n , i.e., 

p w (x) - ft' 2 6{En - ^ (X)) (24) 

Investigations of the semiclassical corrections to this limit have shown that the unstable pe- 
riodic orbits of the chaotic classical system strongly influence the quantum eigenfunctions, 
resulting in accentuations of probability in the vicinity of periodic orbits f§J2^|- Semiclassi- 
cal corrections have also been studied for the stationary Liouville eigenfunctions of regular 



systems [fl5 |. However, no results have been obtained for the correspondence of the nonsta- 
tionary Liouville eigenfunctions or eigenspectrum for the chaotic case. The reason for this 
is made clear in the following paper, where we show that the p^ m (x), n ^ m, for chaotic 
systems do not have classical limits. Rather, one must deal with the classical and quantum 
spectral projection operators discussed below. 

III. CLASSICAL LIOUVILLE EIGENFUNCTIONS AND EIGENSPECTRUM 



The classical Liouville eigenvalue problem [Eq. flTBp] requires that one find a complete 
set of orthogonal Liouville eigenfunctions pa,q( x ) with eigenvalues A, i.e., they must satisfy 
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jdxp*, Q ,(x)p A , Q (x) = (25) 

for A ^ A' or a ^ a', and 

I>A, a ( x o)PA,a(x) = 5(x-x ). (26) 

[Note that for notational simplicity, classical eigenvalues in the remainder of the paper are 
denoted A, rather than A c ]. For the case of s degrees of freedom and r constants of the 
motion we must find simultaneous eigenstates of the complete set of 2r operators [Eq. (0)]. 

If a distribution /(x) is an eigenstate of a real multiplicative operator A(x) with eigen- 
value A', then L4(x) - A']/(x) = 0. That is, /(x) = for all x not satisfying A(x) = A' 
so that all nontrivial solutions /(x) must be proportional to 5(A(x) — A'). That is, all 
distributions which are eigenstates of the multiplicative operators H(x), Ki(x), . . . , if r _i(x) 
must be proportional to 

6{E-H{x))H£6(K j (x)-K' j ). (27) 

Little more can be said about the general case of arbitrary r. We therefore focus on the 
two interesting special cases: r = s, corresponding to integrable systems, and r = 1 which 
includes the case of chaotic motion. 



A. Integrable Systems 

For r = s we can introduce action variables I and conjugate angle variables such that 
H(x) = H(T), and Kj{x) = Kj(T) for j = 1, . . . , s — 1. It follows that there exists an V such 
that 

5{E - F(x))n^- 1 1 5( J ftT i (x) - K'j) oc (J(I(x) - I') (28) 

and that the classical Liouville eigendistributions p\ tCe (x) must be of the form 5(I(x) — 
I')F(0(x)). Since L c = ■ -j^g the Liouville eigenequation is 

- za,(I') • ±F(0) = XF(0) (29) 
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with solutions F(6) = e ia ® , A = a ■ w(I'). If the system is periodic in the angle variables 
then a = k £ Z s , i.e., a is a vector of s integers. Thus the classical Liouville eigenfunctions 
for an integrable system are given by 



with eigenvalue 



Pi', k (x) = " IK 1 ^, (30) 



A r , k = k-u,(I'). (31) 



Note that these eigenfunctions are also eigenfunctions of the differential operators 
i{Kj(x),-} = —i 9K gj^ ■ with corresponding eigenvalues k • dK Qy - ■ If the Kj(I) = Ij 
then the eigenvalues are simply kj. This set of eigenfunctions is complete and orthogonal, 
i.e., 

£ / dV p* llt k (x')pr, k (x) = 5(x - x'), (32) 
k J 

and 

J dx pi„ ik ,(x)p Fik (x) = 5k,k'5(I' - I")- (33) 

We comment briefly on the character of the Liouville eigenfunctions for integrable sys- 
tems; a detailed discussion and related computations are provided elsewhere |fLT| . The 
spectrum of an integrable system on a given torus I' = I(x) is discrete and is given by the 
set of frequencies {k • w(I')|k £ Z s }. On the energy surface E = H(x), the spectrum of 
an integrable system is the union of all sets {k • u>(I')|k £ Z s } from tori I' = I(x) with 
H(I') = E. Typically, w(I') varies with I', so that the spectrum of an integrable system 



on an energy surface is continuous. In our terminology |2T] eigenfunctions with frequencies 
Ai', k = k- cj(I') with k ^ belong to the continuous spectrum while those with k = belong 
to the point spectrum. Examination of Eq. (|30| ) reveals that the eigenfunctions pi',k(I, 0) 
are identically zero off the torus I(x) = I'. Eigendistributions with k = are stationary and 
uniform on the torus, and as a set they project out the longtime limit of any given initial 

10 



probability distribution. For k 7^ the character of the eigendistribution is determined by 
the underlying orbit structure of the torus I(x) = I'. If the frequencies uj 2 , . . . ,u s of the 
torus are commensurate then the orbits of the torus I(x) = I' are all periodic p^j . The 
eigendistributions associated with commensurate tori are nonuniform and stationary when 
k • uj(V) = and k ^ 0, and nonuniform and nonstationary when k • cj(I') 7^ 0. 

The eigendistributions associated with incommensurate tori are nonuniform and non- 
stationary for k 7^ 0. Eigendistributions for integrable systems are supported only over a 
given torus I' = I(x), but every orbit of the torus contributes to the construction of the 



eigenfunction [23 



B. Chaotic Systems 

Now consider the case where r = 1. In addition to being chaotic we assume that the sys- 
tem is hyperbolic, i.e., that it has a countable number of isolated unstable periodic orbits on 
each surface of constant energy. Thus, our chaotic system has positive Kolmogorov entropy, 
exhibits sensitivity to initial conditions (i.e., possesses 2s — 2 nonvanishing Liapunov expo- 
nents), and is mixing, and hence is weak mixing and ergodic ||12|| . Such familiar properties 



of chaotic systems (i.e., ergodicity, weak mixing, mixing etc.) may be related to the spectral 



properties of the Liouville operator [|L2 |. For example, if A = is a nondegenerate point 
eigenvalue of the Liouville operator L c on the energy surface then the system is ergodic. 
If A = is the only point eigenvalue of L c on the energy surface then the system is weak 
mixing. If the rest of the spectrum is continuous then the system is mixing. In this section 
we determine the form of the generalized eigenfunctions which correspond to the point and 



continuous Liouville eigenspectrum ||T2| , [21| for chaotic systems. 

In accord with Eq. ( |TTD we seek eigenfunctions of the complete set L c and H(x). As 
eigenstates of the energy multiplication operator H(x) they must be proportional to 5(E — 
H(x.)) where E is the corresponding energy eigenvalue. Consider then the function r(x) 
which is conjugate to the constant of the motion H(x), i.e., {r(x),if(x)} = 1. Then the 
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distribution S(E — H(x.))e tXT ^ is an eigenfunction of the classical Liouville operator with 
eigenvalue A since L c = i{H(x), ■} = —i-§p- Note that e~ lic *r(x) = r(X(x, —t)) = r(x) — t, 
where X(x, —t) is the phase space point from which x evolves over the time t. Thus, for a 
given choice of the origin of time along each orbit of the system, we can construct a function 
t(x) which specifies the time of an arbitrary point x along the orbit, since trajectories do 
not intersect one another in phase space. However, the distributions 5(E — H (x))e lAr ( x ) ^EA 



do not form a complete set of Liouville eigenf unctions since they do not account for the 
infinite degeneracy of the classical Liouville spectrum for chaotic systems flT2"| . 

To properly define eigenfunctions we must consider the manner in which points in phase 
space are labeled. Specifically we introduce a new system of variables: H(x) which is the 
energy of a point in phase space, t(x) a function which assigns a time to every point in 
phase space, and 2s — 2 stationary variables rj(x) which label the trajectories on the energy 
surface. The functions r and t] are continuous only along an orbit: r increments smoothly 
while r) remains constant along an orbit. The 77 variables are chosen to have the property 
that 

5(x - x) = 6{H(xo) - (x))5(r(xo) - r(x))5(77(xo) - ij(x)). (34) 

As one moves off a trajectory the variables r and 77 must take on all possible values infinitely 
often since (a) there are an infinite number of chaotic trajectories, (b) every chaotic trajec- 
tory comes arbitrarily close to every phase space point on every other trajectory, and (c) 
chaotic trajectories return arbitrarily close to themselves infinitely often. We now introduce 
a complete set of square integrable orthonormal functions Xt(j))i ^ G Z, i.e., 

£xj(l(xb))jfcfa(x)) = *(»l(xo) - r,(x)) (35) 

e 

and 

J dr i xt-ivMii) = Si>4- ( 36 ) 

Classical Liouville eigenfunctions for chaotic systems p E A (x) are then defined as: 
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PU X ) = 4r^ E ~ H^))e^ X ^)). (37) 

V Z7T 

Here the integer I labels the infinite but countable degeneracy of the Liouville eigenvalue 
A for the continuous part of the spectrum [12] at energy E. The distributions p E A (x) are 



supported over the entire energy surface E = Hfe) and so are in some sense constructed 
from every orbit of the energy surface. It can be readily verified that these eigenfunctions 
are both orthogonal and complete (see Appendix A): 

J rfx p^ )A ,(x)4 A (x) = 8 U ,8{E - E')8{\ - A'), 

poo poo 

Y, dE C /Apf iA (x )p^ A (x)= ( 5(x-xo). (38) 

^ JO J— oo 

Here we note that the introduction of the complete set of xe is necessary to achieve the 
correct <5(x — xq) term in Eq. and that these eigenfunctions are quite different from 



those proposed elsewhere which we have shown to be incorrect [17|. 

The eigenfunctions in Eq. fl3"T|) have a structure which at first sight appears to violate 
several well known properties of chaotic systems. For example, the existence of stationary 
square integrable functions x^(t/(x)), which are not functions of the Hamiltonian, appears to 
contradict the fact that for ergodic systems the only stationary square integrable distribu- 
tions are functions of the Hamiltonian. However, the xeivt*-)) distributions are not square 
integrable on the energy surface. That is, 

J dHdrdrj 8(E - H) |^(r/(x))| 2 = J dr = oo. (39) 

In addition, it might appear objectionable that x^(t/(x)) is a global constant of the motion 
independent of the Hamiltonian. However, this is not the case since the set of points x 
satisfying xei'ni'*)) = x\ where x' is some constant, at most contains points from a countable 
number of trajectories and so cannot divide the phase space in any meaningful way. Thus 
the functions xeivi*)) are n °t true constants of the motion. 

The relation between these eigenfunctions and the Liouville spectrum on the energy 
surface E = H(x) is not immediately clear. Examination of the spectral decomposition 
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in Eq. fl38|) reveals no separation of point and continuous spectrum, although the Liouville 
eigenfunctions have the correct degeneracy for the continuous spectrum. A proper treatment 
of the separation of point and continuous spectrum is given in Sec. |V B. 



IV. CORRESPONDENCE: INTEGRABLE SYSTEMS 

Having established the essential features of the classical eigenstate picture we now ex- 
amine correspondence. We assume that r' = r, i.e., that the complete set of commuting 
quantum observables has as many members as there are constants of the classical motion. 

Given the complete set of commuting observables H, K\, . . . , K r -i for r < s, the complete 
set of 2r commuting superoperators is given by [Eq. (|5|)]: 

L, n, l[k lt •], ifc •]+, . . . , •], •]+■ (40) 

We denote the simultaneous eigenstates of the r observables by |n) where n = (ni, . . . , n r ) 
are the quantum numbers, and the corresponding eigenvalues are E a ,Ki(n), . . . ,K r _i(n). 
We further denote the Wigner-Weyl representation of the quantum observables H , Kj by 
H(x), and Kf(x), for j = 1, . . . , r - 1. 

The Liouville eigenstates m must be an eigenstate of a corresponding set of 2r super- 
operators in the Wigner-Weyl representation which are 

L(x), W(x), |A7(x) sin(^), Kj{x) cos(^), . . . 

. . . , jK_,(x) sin(^), K^(x) cos(^), (41) 

with L(x) = (2i/h)H(x) sm(ha/2) and 7Y(x) = H(x) cos(ha/2). By comparison the classical 
complete set is 

L c (x), H(x), z{#i(x), •}, A^(x), . . . , i{K r ^(x), •}, A^i(x)- (42) 

We assume that lim^o KJ{x) = Kj(x). Since \im h ^ cos(ha/2) = 1 and 
lim/^o fffj sin(^o"/2) = «<r, we note that formally the quantum operator set [Eq. ([41])] 
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becomes the classical operator set [Eq. fl42|) l in the classical limit. The first eigenequation 
is the time independent von Neumann (quantum Liouville) equation 

£(*Km( X ) = A n,m^, m (x), (43) 

where £(x) is the quantum Liouville operator in the Wigner-Weyl representation and A nm = 
(E n — E m )/h. Second, p„ m must satisfy the eigenequation for the Hermitian energy operator 
ft = ] + , i.e., the anti-commutator equation 

|n)(m|] + r = ff(x) cos(^/2)p- m (x) = ^^p-, m (x), (44) 

or 

H(x)p- m (x) =J E; n , m p- m (x), (45) 

where E n m = (E n + E m )/2. Similarly, for j = 1, ... ,r — 1, eigenfunctions p„ jm (x) satisfy 
the equations 

(~[^, |n)(m|])- = |i^(x) sin(^/2)p- m (x) = ^ (n) ~ ^ (m) p™ m (x), (46) 

and 

|n>(m|] + )« = ^(x)cos(^/2)p- m (x) = ^y^ pl, m (x). (47) 

We can readily explore the classical limit of the 2r quantum eigenequations, Eqs. ([43|) - 
(0) if 

limPn, m (x) -> p„. m (x), (48) 

where Pn m ( x ) is to be determined. However, this assumption is not generally valid. That 
is, for quantum systems whose classical analogs are chaotic (r = r' = 1), as we demonstrate 
in Ref. |1|], individual quantum eigenfunctions p™ m (x), n ^ m, display essential singularities 
as h — > 0, and hence do not have correspondence limits. Rather, the spectral projection 
operators introduced later below [Sect. V] do have correspondence limits. Nonetheless, we 
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can adopt Eq. (|4q) primarily for notational convenience; the arguments which follow can 
be reformulated in terms of quantum and classical spectral projection operators when Eq. 
(||) is invalid. 

We assume that the h — > limit of Ky(x) is K(pc). In particular, if we consider the case 
where the set of 2r quantum operators [Eq. ([HP] gives the classical set [Eq. (pE2])1 in the 
classical limit and that Eq. (|48 ) holds then each of the eigenvalue equations [Eqs. fl43|)-([47)1 



reduces to the set of classical eigenvalue equations for p\ j(X (x) with the same eigenvalues as 
in the quantum eigenequations. If this is the case then the quantum Pn,m( x ) S oes to the 
classical p\ t0t (x), with clearly identifiable eigenvalues and correspondence is established. 
An example of this approach, applied to integrable systems, is provided below. 

A. Integrable Systems: Formal Correspondence 



In this section we first apply this approach to determine the complete set of correspon- 
dence rules for integrable systems. The resulting picture of correspondence for integrable 
systems is exceptionally clear but gives no insight into the approach to the classical limit. 
This is treated in the next subsection. 

In this instance Eq. ( f48|) has been proven by Berry |Tj| and by Jaffe and Brumer [10 



Thus, in accord with the previous section, the eigenfunctions Pn,m( x ) satisfy 2s classical 
eigenequations in the classical limit: 

£ C ( X K,m( X ) = ^n,mPn,m( X )> (49) 

H( X )pl m ( X )=E n , m pl m ( X ), (50) 
-K m (I, 0) = Kj{n) ~ Ki(m) < m (I, 6) = 0, (51) 

tfiOQpU 1 ' 0) = Kj(n) + / j(m V n , m (I, 0) = 0, (52) 
for j — 1, . . . , s — 1. As in Eq. ( p8| ) we conclude that p^ m (x) is proportional to 
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5(£-tf(x))n-^,(x)-^), (53) 

where E = £ n , m , and K] = ^(°)+ g j( m ) . Since H(x) = H(I), and i^-(x) = K^l) for 
j — 1, . . . , s — 1 there exists an I nm such that 

5(£ - iy(x))n^5(ir i (x) - i^') ex 5(1 - I n , m ). (54) 

Hence p^ m (x) oc 5(1 - I„, m )-Fn,m(0)- 

^From Eqs. (^) and ( |5"TD it then follows that 

- ^(In, m ) • ^Q F n,™(°) = Erl ~ n Em F n AO) (55) 

and 

- ^(In, m ) ■ ^.mW = A \ A V n , m (0). (56) 

By WKB |p!o|,15,[l6l i? n ~ ^(I n ) where I n = (n + /3)7i and where /3 are the Maslov indices. 



Similarly we may approximate Kj(n) ~ ffj(l n ). Assuming that 

In,m - In ~ O(Zl) ~ I n , m - I m << I n , m (57) 

it follows that 

£ n - E m H(I a ) - H(l m ) 

h h 

H \^-n,m (In,m In)) H(\ nm (In.m Im) 



h 

~ ^(I n ,m) ■ (n - m) (58) 
by the Taylor expansion of H(I n ) and H(I m ) about I nm . Similarly one can show that 

(l n , m )-(n-m). (59) 



^■(n)-^(m) dKj 



h dl 

Substituting Eqs. @) and © into Eqs. (0) and @ shows that F njm (0) ~ e 4 ( n - m > . 



Thus, up to a normalization factor 

P«(X) - ^ S/2 pi n , ra ,n- m (x) (60) 
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in the limit as h — > 0. 

More specifically, we have to see that p„ m , an eigenfunction of the 2s quantum operators 

L(x), W(x), |a]»(x) sin(^), A-(x) cos(^), . . . 

. . . , fiC-i(x) sin(^), iC-iW cos(^), (61) 

with corresponding eigenvalues 

^(n)- ^(m) K 1 (n)+i^ 1 (m) 

AT B -i(n) - AT a _i(m) ^x(n) + X,_ 1 (m) 
ft ' 2 

approaches a classical limit ^ m (I, 6) = ^ 2 pi n ,m,n-m( x ) which is an eigenfunction of the 
2s classical operators 

L c (x), tf(x), z{^(x), ■}, ^(x), . . . , z{if s _i(x), •}, ^_!(x), (63) 

with corresponding eigenvalues 

dK 1/T , K 1 (n) + K 1 (m) 



^In : m,n-m, -E-n.m, (n — m) • -^p(I n ,m /• 



91 v w ' 2 

dK a - lf AT fl _i(n) + if s _i(m) 
...,(n-m) .— ^— (I nim ), . (64) 

Using a completely different argument which employed primitive WKB wavefunctions 
Berry showed ]15| that the stationary Liouville eigenfunctions have the correspondence limit 



Pn,„(x) - ^ /2 Pi n ,o(x), (65) 

where I n = [n + (3]h. Since I nn = I n both results are in agreement. For the case of 
nonstationary Liouville eigenfunctions Jaffe and Brumer JK| proved the general formula, 
Eq. flBTf ), where they surmised that I nim = (I n + I m )/2. However, I nim is obtained correctly 
by solving the s equations 

(£n + £m)/2 = #(I n , m ), (66) 

and 



(^(n) + ^(m))/2 = K,(I n , m ) (67) 



for j = 1, . . . , s — 1. The surmise of Jaffe and Brumer |T0[ is an approximation which is good 
to first order in h. To see this note that for cj(0) ^Owe can expand H(T) ~ ^(0) -I + 0(I 2 ). 
Equation (|66"D then reduces to u;(0) • [I n ,m — (In + Im)/2] + 0{h 2 ) ~ for which I nim ~ 
(I n + I m )/2 is clearly a solution to first order in h. If dKj (0) jd\ ^ for all j then it can 
be shown in a similar fashion that I njm ~ (I n + I m )/2 is a first order solution of Eqs. 



Thus, the approximation [[K]] I n jm ~ (I n + I m )/2 is good to first order in Planck's constant. 

B. Integrable Systems: Approach to the Limit 

The method described above assumes the validity of Eq. ( fjE8"| ) and fails to provide insights 
into the behavior of the eigendistributions as h — > 0. Here we extend and reexamine our pre- 
vious proof |Tt| of Eq. (^) for integrable systems with particular interest in the mechanism 



by which essential singularities are avoided as h — > 0. The following argument also deter- 
mines the correspondence limits of the stationary and nonstationary Liouville eigenfunctions 
in a unified manner. 

In the h — > limit the eigenfunctions of the Hamiltonian, (q|n), are well approximated 
by linear combinations of primitive WKB wavefunctions, i.e., 

(q|n) ~ ^(27T)- g /Vet 92 ^ ( ^ In) r /2 e^ (q ' In)/?i e-^ 7r/2 (68) 

where i?j(q, I) = Jq o pj(q',I) • dq' are the actions corresponding to the different branches 
of the multivalued momentum functions pj(q, I), and Uj are the Maslov indices. Thus, it 
follows that 

Pn, m (x) = hr^\ dv C W* <q _ v/2|n)(m|q + v/2) ~K(x), (69) 

where 



3,k 



C&(*) = (2nh)-^ J d^\det d2R ^-J^ In) det ^^J^ ^ | 1/2 ■ 

e ip-v/h e iRj(q-v/2,I n )/h e -iR k (ci-+v/2,l m )/h e i(uj-v k )TT/2 ^Q'j 
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and where, in principle, we must retain all contributions to the double sum from the different 
sheets of the phases. However, one of us has argued |TIJ that the j ^ k terms should 
make a negligible contribution to p„ m (x) due to the rapidly oscillatory character of their 
integrands. Before examining the j — k terms which yield the classical limit as h — > 0, we 
briefly demonstrate that the j ^ k terms are indeed small if some averaging over phase space 
(essentially an energy average) is performed. (We would prefer to carry out an average over 
wavefunctions at different energies, as in the chaotic case 0, but it is unclear as to how to 
do this. We expect the phase space average to achieve a similar result). 

Making the change of variable v — > hv Eq. fl70|) can be written in the form 



(71) 



&> k fx) - h s ' 2 f ivldct 92jRj(q ~ W2 ' In) dct g2jRfc(q + W2 ' Im) 1 1/2 ■ 
^W-n Jdv|det ^ det d ^ Ql | 

2-iripv iRj(q-hv/2,I n )/h -iR k (q+hv/2,I m )/h i(fj -v k )ir /2 



Note the presence of essential singularities in each of the factors e lR ^ hv / 2 M/ h and 
e -iR(q+hv/2,i m )/h _ E X p an ding the actions R in powers of h, we have 

h h 2 dr) (a I ) 

R 3 ( q - hv/2, I n ) ~ Rjid, I n ) - -p,(q, I n ) • v + -v ■ ° PA ^ n> ■ v, (72) 

and 

R k (q + hv/2, I m ) ~ R k (q, I m ) + ^p fe (q, I m ) • v + • . v . (73) 

We also assume that the determinant factors in the integrand are slowly varying, i.e., 

ldct d 2 Rj(q. - hv/2, I n ) d 2 R k (q + hv/2, I m ) 1/2 ^ ^ ^-(q, I n ) ^ ^^(q, I m ) 1/2 
<9q<91 <9q<91 <9q<91 <9q<91 

(74) 

Here the tilde denotes the transpose (a column vector) of the row vector v. Substituting 
these expressions into Eq. ([n]) and performing the integral over v shows that Eq. (|71~D 
becomes 

Q fc m (x) = (20r)'/ 2 |det d6t d2R d^di m) ' V2 ^ ~ gP^0 r i/2 . 
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where p^k = (pj(q, I n ) + Pfc(q, I m ))/2. The contributions C^m(x) are therefore highly 
oscillatory functions of the phase space variables. 

We now smooth Eq. (f73|) over intervals of length Ap about each momentum variable pi, 
I — 1, ... s. In the stationary phase approximation (for h — *■ 0) 

• f ^ A, r^ dp 3 (q,I n ) _ %(q,I m ) ,_ 1/2 _ 



J~Ap/2 aq aq 



exp{-2?(p + u - p j>fe ) • [ — ^ ] • (p + u - Pj, k )/h} 

-> (7r^) s / 2 fi Ap (p- Pj - fc ), (76) 



where Oa p (p) = 1/ (Ap) s if p^ G [pi — Ap/2, pz + Ap/2] for / = 1, . . . , s, and is zero otherwise. 
It follows that smoothing Eq. ([75]) over a classically small interval in momentum gives 



C^ m (x) ~ h'' 2 \fa ¥yj^ - pj, 



(77) 



Further smoothing C^ k m (~x) over intervals of length Aq about each coordinate variable 
qi, I = 1, . . . s, one can show that the C^ k m (x) contribution to /}JJ m (x) is proportional to 

W l=1 smcU-p l k )Aq/2h] (78) 



where p^- is the I th component of p 3 -(q, I). Thus, for j ^ k Eq. (|78| ) is 0(h s ), while for 
j = k Eq. (fTTf ) is 0(1). Thus, we have shown that smoothing over phase space effectively 
eliminates contributions from the j ^ k terms. 

Focusing now on the j = k terms, and dropping the subscript j, it follows that the 
diagonal contributions take the form 

C£L(x) ~ ^Idet *^^ det 9 ^ ( q q 9I Im) 1 V2 e iW q ,i n )-^ (q ,i m ))/^ (p _ 1 [p(q> In) + p(q Im)])> 

(79) 

where we have replaced ^A p (p— §[p(q, I n )+p(q, Im)]) by the delta function S(p— |[p(q, I n ) + 
p(q, I m )]) since we are free to choose Ap small. Assume now that there is an action I nm 
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such that p(q, I n ,m) ~ |[p(q? In) + p(q? I m )] ( see Ref. 0]), and we also assume that I n , I m , 
and I nim satisfy relation (|57D , then it readily follows that 

R(q, I n ) - R(q, Im) R(<1, In,m ~ (I n ,m ~ In)) ~ R(<1, I„,m ~ (In,m ~ Im)) 



dR(q, I n , m ) 



n — m 



(80) 



dl 

by Taylor expansion of i?(q, I n ) and i2(q, I m ) about I n)m . Further since 9R ^'^ = 0, we 
obtain 

Q,Ux) ~ ^[detgg^det ^^ ^(p - p(q,I nim )) e ^)- W. (81) 

Finally, summing over the different sheets j [Eq. flB9p], assuming that for small h 

| det ^M det ^|W|V, „ |d**^£*>|, (82) 
aqol aqol oqal 

and changing variables in the delta function we have 

Pn,m(x) ~ ^ /2 5(I(x) - I D , m )e««-™>°M, (83) 

which is the correct correspondence limit. In particular, p^m becomes the classical Liouville 
eigenfunction [Eq. (|30|) ] with I = I n m and k = (n — m) in the classical limit. Two things 
are noteworthy. First, an average over ApAq was necessary to eliminate highly oscillatory 
terms in the j ^ k contribution to Eq. (p9|) . Second, under these circumstances the essential 
singularities in the original expression cancel. 

V. CLASSICAL LIOUVILLE SPECTRAL DECOMPOSITION 

Attempts to follow a similar correspondence approach for chaotic systems fails. That 
is, Eq. (HD is not true insofar as individual eigendistributions p™ m (x) do not have proper 
classical limits. Examination of our motivating equations fll2|) and (|20"D suggests that a one 
to one correspondence of quantum to classical eigenfunctions is not necessary to establish 
quantum-classical correspondence as h — > 0. Rather, we require a relation between the 
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quantities in brackets in Eq. (|12|) and fl2~0|) as h — > 0. In this section this idea is quantified, 
for classical systems, by introducing the classical spectral projectors which we subsequently 
relate in the h — *■ limit, to the analogous quantum projectors. In addition, we provide 
a general approach to computing the classical projectors for arbitrary systems. 

The classical projectors will be seen to offer a better method of describing chaotic systems 
than do the classical Liouville eigenfunctions introduced above [Eq. ([371)1. Specifically, they 
are not based on the arbitrary set xt, they are completely defined as opposed to the p E A (x) 
which are defined to within an overall phase, and they offer a means of separating the 
spectrum into continuous and discrete components |[21|| . 

We also extend the approach and separate the singular spectrum of the periodic orbits 
from the rest of the continuous spectrum. The resulting periodic orbit spectral projection 
operators are used in the companion paper on correspondence in chaotic systems 0. 

To develop an energy /frequency spectral decomposition we must isolate the fixed fre- 
quency and energy components of an arbitrary time evolving phase space probability density 
pipe, t) for a conservative classical Hamiltonian system of s degrees of freedom. Phase space 
probability densities p(x, t) = p(X(x, —t), 0) evolve in classical mechanics in accord with the 
classical Liouville equation, Eq. flT5p. Consider that 

p(x, t) = J rfx p(x , 0) 5(x - X(x, -t)) (84) 

where 5(x — X(x, — £)) is the kernel of the Liouville propagator, which can be rewritten as 

/OO P 
dt'S(t - t') / rfx p(x , 0) 5(x - X(x, -t')). (85) 

Making the replacement 5(t — i') = dX e-^it-t') we Q^^ain 

/oo 1 roc r 

dt > / dXe-^ 1 - 1 ) / dx p(x , 0) 5(x - X(x, -t')). (86) 
-oo Z7l J— oo J 

Finally inserting a closure relation for the energy 1 = / °° dE 5(E — H(x.)) gives 

roc roc r 

p(x, t) = J o dE J dX e^ xt J rfx p(x , 0) T SiA (x; x ), (87) 
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where 

T s , A (x; x ) = — 6{E - H{*o)) / dt'e iXt 5(x - X(x, -t')). (88) 

Z7T J-oo 

Equation (|87| ) provides an expansion of p(x, t) in terms of the classical spectral projection 
operators Y^a which are expressed, in Eq. (|8|), in terms of the dynamics. Specifically, 
T£j ; a(x; x o) is the Fourier transform of the kernel of the Liouville propagator, restricted to 
the energy shell. Equation fl87|) is a restatement of Eq. (^) which emphasizes the classical 
spectral projection operators as being central to the analysis of the dynamical evolution of 
distributions p(x, t). 

The distribution Y^^xjxo), as a function of x, is an eigenfunction of the Liouville 
operator with eigenvalue A. To see this consider that propagation with respect to the x 
variables gives 

e- iLct T £ , A (x; x ) = —S(E - (x )) / dt'e lXt 5(x - X(x, -i! - t)) 

2,71 J-oo 

= e- iA 'T SiA (x;x ), (89) 
where we have made the change of variable t'+t — > t'. Similarly, with L c (x ) = i 91 ^^ -J--^- 
operating on the x variables, L c Te,a(x;xo) = — ATg^x; x ). [Here J = 



is the 



\< / 



2s x 2s dimensional symplectic matrix [25[]]. Thus, T^ A (x;x ) is an eigenfunction of L c 
with eigenvalue A (—A) in the x (x ) variables. These properties are consequences of the 
fact that T^ ) a(x;x ) = T^_a(xo;x), which is in turn a consequence of the fact that L c is 
unitarily equivalent to — L c . 

The decomposition in Eq. (|87| ) holds for all distributions and all phase space points x, 
including periodic orbits. To show this, suppose that x lies on a periodic orbit. Then by 
Eq. 



T s , A (x; x ) = —5{E - (x )) / dt'e iXt <5(x - X(x, -t')) 

Z7T J-oo 

= 6{E - H(xo)) E ^' Ar / dt'e iXt 5(x - X(x, -t')) 

2lT , ' J-t/2 



J=-oo 
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8{E - H(xo)) £ 5 ( X ~ 1*3 rt- / dt'e 2 * ijt '/ T 6(xo - X(x, -t')) (90) 



j=~oo 

where r is the fundamental period of the orbit. Thus we have 

1 rl 2 



T £iA (x; x 

j=-oo 



8{E - F(xo)) £ *(A - 27rj/r)- / dti*™* ' T 8{x* - X(x, -0) (91) 

on periodic orbits. Now inserting this expression into the right hand side of Eq. ([87]) and 
performing the integrals over A and E gives 

1 r/ 2 



p(x, t) = / rfx p(x , 0) £ - / dH^^S^ - X(x, -0) 

/• 00 /-t/2 

= yrfx oP (x ,0) £ y^ 2 dt , 5(t / -t-jr)5(x -X(x,-t')) 



j=-oo 



= I rfx p(x , 0) 5(x - X(x, -*)). (92) 

Performing the integrals over <ix gives the probability density on the point x of the periodic 
orbit, as required. Thus, the expansion fl8T|) and its associated closure relation hold for points 
on periodic orbits. 

Since Eq. fl57p holds for all p and all x it implies the existence of the closure relation 

/*oo /*oo 

dE d\ T £ . A (x; xo) = 5(x - x). (93) 

JO J -co 

Two unsatisfactory issues require more careful attention: First, there is no clear separation 



of the point spectrum from the continuous spectrum [|2T|]. Second, a comparison of Eq. 
with Eq. (p0|) suggests that T S A (x, x ) should relate to the Liouville eigenf unctions p E x as 
a product of an eigenfunction in x times an eigenfunction in x. These issues are addressed 
below where we treat the integrable and chaotic cases separately. 



A. Integrable Systems 

To demonstrate the computation of T consider integrable systems. Here we change to 
action angle variables x = (I, 6) and x = (Iq, ) so that Te,\ assumes the form 
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1 



T £ . A (x; x ) = —S(E - H(I )) / dt'e m 5(I - I)6((8 -6 + w(I)t') mod 2tt) 

= - tf(Io)) T ^e^'5(Io - 1)-^- E '^W, 

Z7r j-oo [Zir) kgZs 



where we have used the identity |26| 



E *(*o - + "(I)*' - 27TJ) = j7T-r s E 
jGZ s v Z7 v kez s 



3 -ik-[0 Q -0+u;(i)i 



In terms of the Liouville eigendistributions pi> >n (l, 0) [Eq. (^)] Eq. (|94j) becomes 

T £iA (x; x ) = - H(I )) E / ^I'5(A - n ■ u>(I'))p* v n (I , O )pi>(I, 0). 

2vr n J 

Inserting Eq. (|96|) into Eq. (|87|) gives the expansion 

p(x,t) = E / di'e-™-"^ ( dl d0 P (i ,e ,o) k n (i o ,0 o )pi>(i,0 



(94) 



(95) 



(96) 



(97) 



This result suggests a natural definition of spectral projection operators for integrable sys- 
tems, i.e., 



^r,n(x;x ) = pi, n (x )pr in (x), 



(98) 



which takes the form of a product of Liouville eigenfunctions. These spectral projection 
operators satisfy the following symmetry, closure and orthonormality relations: 



^r,n(x;x ) = 3 ; i'- n (x ;x) 



(99) 



and 



J dl' y Fi o(x; x ) + E / dl' ^r, n (x; x ) = 5(x - x ), 



J rfx ^(x;^) ^i>(x;xo) = 5^8(1" - I / )^, n (x[ ) ; x ). 



(100) 



;ioi) 



Further, the decomposition in Eq. ( |100| ) clearly displays the separation of point and contin- 
uous spectrum. 
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B. Chaotic Systems 



The time evolution of a phase space distribution for a chaotic system which is hyperbolic 
can be decomposed in the following fashion (see Appendix B) 

p(x,£) = J dE J cfx p(x ,0) T#(x;x ) + J dE j-dX e~ lXt J dx. p(x ,0) T SjA (x;x ), 



(102) 



where Te^(x;x ) is given by Eq. (|88|) and 



T,fx- x 1 - 5{E - H ^^ E ~ Q 03) 

Tb(Xi Xo) " fdx>5(E-HW) ' ( } 

Here -fdX denotes the integral over A with the point spectrum at A = removed, i.e., 

-fdX = f^co dXPE, where PeP = [p — (p)e], and 

Jdx5(E-H(x)) p(x,0) 

(p)s " ;^(E-i/(xO) (104) 

is the micro canonical average of p. The operator Pe projects distributions onto the contin- 
uous spectrum part of the Hilbert space. Properties such as the system is ergodic (A = is 
a nondegenerate point eigenvalue with a corresponding eigenfunction which is uniform over 
the energy surface), weak mixing (no point eigenvalues other than A = 0), and positive Kol- 
mogorov entropy (i.e., the remainder of the spectrum is continuous) have been incorporated 
in Eq. ( |102| ) in an obvious way. Since the component of p(x) with energy E is given by 



J dx p(x , 0)T £ (x; x ) = (p) E S(E - H(x)) (105) 

we see that projects p(x) onto the point spectrum part of the Hilbert space. Similarly, 
Te,x projects an arbitrary distribution onto the continuous spectrum part of the Hilbert 
space associated with a frequency A and energy E. We further note that this set of projection 
operators is complete 

pea pco /• 

/ dEr E (x;x )+ d£|rfAT £A (x;Xfl)=i(x-x ), (106) 
Jo Jo J 

and orthogonal: 
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J dx T^(x; x' ) T E (x; x ) = 5(E - E')T E (x.' ; x ), (107) 

and 

J rfx T^ A ,(x;x' ) T £iA (x;x ) = 6(E - E')5(X - A')T £ , A (x[ ) ; x ). (108) 

[These relations are proven in Appendix B.] The distributions T E and T^,a are thus the 
kernels of orthogonal stationary and nonstationary spectral projection operators for the 
classical Liouville spectrum. We have thus achieved, in this approach, the separation of 
point and continuous spectrum. It remains to examine the relationship of Y E> \ to the 
Liouville eigenfunctions. 

To do so insert Eq. fl34|) into the definition of T E ,x [Eq.fl5B|)1 to obtain 

T B , A (x; xq) = - H(x ))8{E - H( X ))e- iX ^-^5(ri(x ) - ry(x)). (109) 

Ztt 



Inserting Eq. (|53| ) in Eq. ( |109| ) and rewriting in terms of the definition of the chaotic 
Liouville eigenfunctions [Eq. (P7|)] gives 

T £ , A (x;x ) = £p£, A (x )p^ A (x). (110) 

i 

We also define stationary Liouville eigenfunctions, Pb(x), via 
such that 

T£;(x;xo) = p^(x )p£;(x). (112) 

The distributions pe(x) and p|; A (x) are zero off the energy shell E = H(x). The eigendis- 
tributions p E (x) are stationary and uniform over the energy shell and belong to the point 
spectrum. Distributions p^ A (x) are nonuniform and stationary for A = 0, and nonuniform 
and nonstationary for A ^ 0. The distributions p|; A (x) belong to the continuous spectrum. 
Finally, note that p.e(x) and pg A (x) are supported over the entire (2s — l)-dimensional en- 
ergy surface E = -ff"(x), in contrast to the distributions pi',k( x ) for integrable systems which 
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are supported only on the s- dimensional torus I' = I(x). In other words, every orbit of the 
energy surface E = H(x) contributes to the construction of a chaotic eigenfunction. 
Comparing 

T B (x; x ) = p^(x )p£;(x) (113) 

with 

r £ , A (x;x ) = Ep^XoVi^x) (114) 

i 

we see that the stationary spectral projection operators are simple products of Liouville 
eigenfunctions while T#,a are composed of a sum of products of Liouville eigenfunctions. 
This structure, i.e., a sum over products, has important implications for the correspondence 
of the nonstationary chaotic quantum Liouville eigenfunctions [|IJ . Hence we shall prove here 
that T SjA (x; x ) cannot be rewritten as a simple product of some new, as yet undiscovered, 
eigenfunctions. 

Suppose that we can write 

T s , A (x; xo) = Q^ >A (xo)Q EiA (x) (115) 

for some unknown distributions $7e )A (x). Then, in order to satisfy the orthogonality relation 
[Eq. (|108|) 1 we must have 



J rfx n^ iV (x)fi B ,A(x) = S(E - E')S(\ - A'). (116) 

Equations ( |115| ) and (|116|) along with the definition of Te,\ then imply that 

8{E - E')8{\ - A')^, A (x) = — 6{E - F(x)) / dt'e iXt fi^y (X(x, -t')) 

= 5{E - H{x))6{\ - A')^,A'(x). (117) 

Integrating over A' and E' gives 

^ >A (x) = 5{E - #(x)) J dE' ^, A (x), (118) 
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which implies that Qe,\(x) oc 5{E — H(x)). Note that Te,a(x;xo) = unless xo and x lie 
on the same trajectory. Thus, if Qe,\(*-) = S(E — iif (x))iV\(x) then iV^(x )A^ A (x) = unless 
xo and x lie on the same trajectory. However, if xo and x lie on the same trajectory then 
iJ(x) = i?(x ) and so 

^ jA (x )fi E , A (x) = 8{E - H(x ))N* x (x )5(E - (x))JV A (x) 

= 6(E-H(jc)) 2 N*(x q )N x (x) (119) 

so that Q* E x (~X())Qe,\{x) would not even be integrable in a delta function sense. Thus Te,\ 
cannot be written as a simple product of Liouville eigenfunctions. 



1. Periodic Orbits 

Periodic orbits play an important role in modern theories of quantization in chaotic 
systems. Given that Eqs. ( 87|) and ( |102| ) hold for points on periodic orbits, and given that 



hyperbolic systems have a countable number of isolated periodic orbits for any given energy, 
it is also possible to separate the singular spectrum associated with the periodic orbits 
from the rest of the continuous spectrum. To do this we define periodic orbit distributions 
(extracted from Eq. (|9lD): 

T* (x; x ) = 5(E - # (x ))- dt'e^'^d^ - X(x, -t')) (120) 

Tk J-T k /2 

where x fc is a point on the k th periodic orbit at energy E. [Here again is the fundamental 
period of orbit k] . Note that the distributions T E have the property that 
/cix^xp T| [0 (x;xo)F(xq)G(x) fc 

Sdx ,6(E-H(X)) = WbWw (121) 

where F and G are any two observables, (F) e denotes a microcanonical average at energy E, 
and (G) Ej = — Jl k J k 2 / 2 dt'e 2mjt '^ Tk G(X(x fc , t')) is the j th Fourier component of G on periodic 
orbit k at energy E. The total closure relation now takes the form 

/ dE T£;(x;xo) + dE S rj^ ),r] k ^E,j(^ x o) + / dE -J-dX T £iA (x;x ) 

= 5(x-Xq), (122) 
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where V = {fl\r) labels a periodic orbit} and the rj variables are those which we introduced 
earlier, and where -fdX now denotes the integral over A with the point and singular spectrum 
removed, i.e., -fdX = J dX SePe, where 

S E p(x) = p(x) -EE <%x),*7 fc (p)ej (123) 

and so 

S E P E p(x) = P(x) - (p) E -EE *K*),i?> " (p>b]>Ij- (124) 

The distributions -(x;xq) are asymmetric in the x, x variables, and so a Kronecker 
delta 5r7(xo),T7j, is introduced as a factor in the singular spectrum term of Eq. (|122|) in order 



to restore the proper symmetry with respect to x and x . The prime over the sum on j 
denotes that the point spectrum has been removed, i.e., Ylj — J2j Pe- 

One last aspect of definition ( [L20|) should be emphasized. For j = we may write 



T| i0 (x; x ) = S(E - H(xo))- T ' dt' 5(r(p\ q fc ) - r(p, q) + if) ■ 

T k J-T k /2 

5(E - H(x))5(r, k - T7(x)) 
= r.^iE - H(xo))5(E - H(x))6(r) k - i|(x)). (125) 



But local coordinates £fc(x) of the Poincare surface of section [20,27], transverse to periodic 
orbit k, may be introduced such that 

5(77 fe -77(x))=5(£ fc (x)). (126) 

Note that rj k — ry(x) ^ Cfe( x ) since ry(x) are global r independent variables, while £ fe (x) are 
r dependent local variables, i.e., £ k (H, r + Tfe, 77) = M k £ k (H, r, 17) where is the stability 
matrix of orbit k. However, since 

8{£ k {H, r + r k , 77)) = T^-Mk(H, r, 77)) (127) 

and I detAffc | = 1, it follows that 6(£ k (H,T,rj)) is time independent. Thus, the stationary 
periodic orbit spectral projection operators can be written in the form 
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T* f0 (x; xo) = r k l 5{E - H(^))S(E - (x))5(£ fe (x)). (128) 

This particular form will be of use in the companion paper |]J on correspondence in chaotic 
systems in which scar contributions are related to these spectral projectors. Indeed, the 
analysis of the classical dynamics of distributions, contained in this section, will prove central 
to the analysis of classical-quantum correspondence in chaotic systems, discussed in the 
following paper 

VI. SUMMARY 

In summary, we have constructed a coherent framework for the study of quantum-classical 
correspondence. We have emphasized the importance of considering, in each mechanics, the 
entire set of 2r eigenequations in the Liouville picture. We have then shown how this ap- 
proach, when combined with previous demonstrations of correspondence for the Liouville 
eigenfunctions, allows a complete understanding of correspondence for integrable systems. 
General methods were then discussed for the construction of Liouville spectral decomposi- 
tions necessary for the study of correspondence in chaotic systems. These methods were 
then employed for the construction of Liouville eigenfunctions for chaotic systems. It was 
also shown that application of these methods to integrable systems yields the usual spec- 
tral decomposition. Finally, we reviewed the arguments for the correspondence of Liouville 
eigenfunctions for quantum systems with integrable classical analogs and showed that the 
primary mechanism of correspondence is the elimination of essential singularities when there 
is an averaging over a small range ApAq. 

In the following paper we show that the correspondence problem for systems with 
chaotic classical analog can be treated by the methods introduced in this paper. 
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APPENDIX A 

Here we will prove the orthogonality and completeness relations ( |38| ) for the chaotic 
Liouville eigenfunctions p% \- To prove ( |38| ) consider that 

/ rfx p& v (x )pi A (x) = / dHdrdr) ±- f°° dre^ x ^ T ■ 

J J Z7T J-oo 

5(E - E')5(H - E) X *Ari)Xe(ri) 
= 5 e/ 5{E' - E)S(X' - A) (129) 

where we have used Eq. fl36[) and the identity 8(x) = ^ Jf^ dy e txy . 
To show completeness note that 

/*oo /*oo ! /*oo /*oo 

dE d\p% x (x>)p f EiX (x) = -Y, dE d\6(E-H(x>)- 

l JO J-oo Z7T J JO J-oo 

S(E - iJ(x) e ^W-^ x ')) X K^(x , ))^(^(x)) 
= 6(H(x) - ^(x'))5(r(x) - r(x'))%(x) - r^x')) 

= 5(x - x') (130) 

where we used Eq. ([$5|) and the identity 8(x) = ^- /f^ dy e lxy to go from the first line to 
the second, and Eq. ( PH ) to go from the second line to the third. 

APPENDIX B 



Here we will prove the decomposition of Eq. (|102|) and relations ( |106|) -( |T0~8l ). To begin 



we will prove relation (|102|) . Let p(x, t) be a solution of the classical Liouville equation. The 
following decomposition can be verified by inspection: 

, \ f°° f , , n J(E - H(xo))5(E - H(x)) 

("OO 

+ / dE p ns {x,t)8{E - H{x)), (131) 
J o 
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where 



p ns (x,t) = p(x,t) - (p) 



(132) 



with (p)e as given in Eq. ( |104j ). Furthermore, we may write 



p ns (x, t) = J dxi p„ s (x , 0)5(x - X(x, -t)) 

/OO P 
dt'5(t' - t) / dx p ns (x , 0)5(x - X(x, -*'))■ (133) 
-OO J 

Substituting i rfA e^^'" 4 ) for 5(t' - i) in Eq. ([133]), inserting 1 = / °° d£ <5(£ - if(x)), 
and then substituting Eq. ( |133| ) into Eq. ( |132j ) and collecting terms, we obtain 

pOO p 

p(x,t) = dE J dx p(x ,0)T jB (x;x ) 

POO P P 

+ J dE jd\e- iXt J dx p(x ,0) T E , A (x;x ), (134) 

where -fdX p = J dX Pe p = J dX [p — (p) e], which proves Eq. (|102|) . Since Eq. (|134|) must 
hold for all t it must also hold for t = 0. We thus have shown that 

poo p poo p p 

p(x, 0) = J dE J dx p(x ,0)T £; (x;xo) + J dE j-dX J rfx p(x ,0) T E>A (x;x ), (135) 

and since this must hold for all p, relation (|106|) must be true. Here of course we have used 
definitions (|103|) and fl88f) . 

The proof of relation (|107|) is much simpler. Inserting definition ( 103 ) into the left-hand 



side of relation (|107|) gives 



J dx T^,(x;x[ ) )T i; (x;xo) = / dx 



5(£' - (x^(£' - (x)) 



[/ dx'5(£-#(x'))] 2 
5(E-H(xo))5(E-H(x)) 

^ ^-g(x^))^-g(x )) 

Simple manipulations and use of definition (|103| ) then give Eq. ( |107| ). 

Relation ( |108| ) can be proven as follows. From definition (|5BD it follows that 

J dx T^, )A ,(x;x[ ) )T £ ; iA (x;xo) = 

— ^S(E - E') I dx r dtdt' 5{E - H{x))e~ iXt e iyt ' ■ 

yZTT j J J — oo 

S(x' -X(x, -t))5(xo - X(x, -f)) 



(136) 



(137) 
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and this can be rewritten in the form 



J dxT^, A ,(x; x' )Te )A (x;xo) 



8{E-E') / dtdt'5{E - H{x ))e- lXt e 



iXtAX't 1 



(2tt) 2 

5(x -X(x' , -(*'-*))) 



138) 



and now changing t' — > t + t' we readily obtain Eq. ( |108| ) 
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